Anomalous hydro dynamic interaction in a quasi-two-dimensional suspension 
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We have studied the correlated Brownian motion of micron-size particles suspended in water and 
confined between two plates. The hydrodynamic interaction between the particles exhibits three 
anomalies, (i) The transverse coupling is negative, i.e., particles exert "anti-drag" on one another 
when moving perpendicular to their connecting line, (ii) The interaction decays with inter-particle 
distance r as faster than in unconfined suspensions but slower than near a single wall, (iii) At 

large distances the pair interaction is independent of concentration within the experimental accuracy. 
The confined suspension thus provides an unusual example of long-range, yet essentially pairwise 
correlations even at high concentration. These effects are shown to arise from the two-dimensional 
dipolar form of the flow induced by single-particle motion. 

PACS numbers: 82.70.Dd, 83.80.Hj, 47.60.+i, 83.50.Ha 



Liquid suspensions containing particles of nanometer- 
to-micron size (colloids) are ubiquitous in industry and 
biology . In various circumstances the colloids are spa- 
tially confined, e.g., in porous media, biological constric- 
tions, nozzles, and microfluidic devices 2]. The dynamic 
behavior of suspensions is affected by flow-mediated ve- 
locity correlations. In unconfined suspensions 1] these 
hydrodynamic interactions are positive (particles drag 
one another in the same direction), long-ranged (decaying 
with inter-particle distance r as 1/r), and involve many- 
body effects (i.e., depend on concentration) P. Consid- 
erable efforts have been devoted recently to studying col- 
loids at the microscopic level, particularly in confinement, 
using digital video microscopy 4] and optical tweezers 
. Particle dynamics near a single wall 0, ^] , between 
two walls 0, in a linear channel 9^, and in a finite con- 
tainer 1^ were investigated. These studies have high- 
lighted flow-mediated effects of the boundaries on parti- 
cle dynamics. Quasi-two-dimensional (q-2D) suspensions 
have been studied recently also by computer simulations 
[in| . In this Letter we demonstrate, using digital video 
microscopy, that the hydrodynamic interactions in a q- 
2D suspension are drastically different from those in less 
confined as well as more confined systems. 

The experimental system consists of an aqueous sus- 
pension of monodisperse silica spheres (diameter 2a = 
1.58 ± 0.04 /im, density 2.2 g/cm^, Duke Scientific), 
tightly confined between two parallel glass plates in a 
sealed thin cell (Fig. H). The inter-plate separation is 
w = 1.76 ± 0.05 /im (llj, i.e., slightly larger than the 
sphere diameter, 2a /w 0.90. Particles thus move 
essentially parallel to the plates, forming a q-2D sus- 
pension. The glass surfaces are coated with chlorine- 
terminated polydimethylsiloxane telomer (Glassclad 6C, 
United Chemical Products) to avoid particle sticking. 
Digital video microscopy and subsequent data analysis 



are used to locate the centers of the N spheres in the field 
of view (area A — 106 x 80 /im^) and then extract time- 
dependent two-dimensional trajectories (time resolution 
0.033 s). Details of the setup and measurement methods 
have been described elsewhere J^J. Measurements were 
made at four packing fractions, = Nno? /A = 0.254, 
0.338, 0.547, 0.619 ± 0.001. From equilibrium studies of 
this system 13] we infer that, for the the current study, 
the particles can be regarded as hard spheres. 




FIG. 1: (a) Optical microscope image of 1.58 /xm-diameter 
silica spheres suspended in water and confined between two 
glass plates at area fraction (j) = 0.338. (b) Schematic view of 
the system and its parameters. 

For a pair of particles in two dimensions, separated by 
a distance r, the coupled Brownian motion can be charac- 
terized by four diffusion coefficients, denoted D^{r) and 
D^{r). The longitudinal and transverse collective coef- 
ficients, and D^, characterize the diffusion of the 
center of mass of the pair along and perpendicular to 
their connecting line. The relative coefficients, and 
D^^ relate to the fluctuations in length and orientation 
of the vector r connecting the pair. These coefficients 
are measured as D^{r) = {[xi{t) ± X2(t)]^)r/(4t) and 
^t(0 = {[yi{t) ±y2{t)]^)r/{4^t), where Xi{t) and yi{t) 
are the displacements of particle i during the time in- 
terval t along and perpendicular to the line connecting 
it to its partner. The average (• • •)r is taken over all 
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observed pairs whose mutual distance falls in the range 
r±0.09 /am. The time interval was kept sufficiently short 
{t < 0.2 s) that r could be assumed constant within the 
±0.09 jam bin size. The pair diffusion coefficients are 
defined above such that, in the absence of correlations 
(e.g., when r oo), all four reduce to the self-diffusivity 
of a single particle, D^. (This is also verified experimen- 
tally.) Thus, the deviations of t fro^^ serve as our 
measure of the pair hydrodynamic interaction. 




2 4 6 

B 



10 12 14 16 18 20 

r(|im) 



0.08 



i 0.07 



0.06 



0.05 



0.04 




12 14 16 18 20 



r (|im) 

FIG. 2: (Color online) Longitudinal (a) and transverse (b) 
pair diffusion coefficients as functions of inter-particle dis- 
tance. Filled and open symbols correspond, respectively, to 
collective and relative coefficients. Symbol shapes and colors 
correspond to different area fractions: (j) = 0.254 (black cir- 
cles), 0.338 (red squares), 0.547 (blue triangles), 0.619 (green 
diamonds). Note the long range of the interaction and the 
sign reversal of the transverse coupling. 



and distances according to 

Al^{p) = {D^^^-D,)/{Doa/w), p = r/w, (1) 

where I^o = k^T / {^Tirjo) is the self-diffusivity of an iso- 
lated, unconfined sphere {k^T being the thermal energy 
and T] the viscosity). This rescaling yields, in the limit 
of small particles and low concentration, parameter-free 
functions describing solely the pair interaction The 
rescaled results are given in Fig.lHl where all long-distance 
measurements collapse onto the same four curves. 
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FIG. 3: (Color online) The same data as in Fig. 2 rescaled 
according to Eq. (QJ. Note the collapse of the long-distance 
data for different area fractions onto the same four curves 
(insets). The solid curves correspond to ±A/p^ with A = 0.31 
[Eq. 0]- The dashed curves in (b) are obtained from the 
rigorous solution in the limit a <^ w ^3 (^^^ text). 



Figure El shows the measured values of t(^) func- 
tions of r at different area fractions As (j) is increased 
the curves shift to lower values (i.e., expectedly de- 
creases with concentration) and, at short distances, be- 
come more structured 14]. We focus in this Letter on the 
long-distance pair interaction. To get a more unified pre- 
sentation in this domain, we rescale diffusion coefficients 



In the longitudinal modes, the correlation has the ex- 
pected positive sign — particles drag one another in the 
direction of motion, thus enhancing the collective mo- 
tion and suppressing the relative motion [Fig.|3Ja)]. The 
transverse interaction, however, has a surprising reverse 
sign — the collective motion is suppressed while the rela- 
tive motion is enhanced [Fig. |3fb)], i.e., particles exert 
"anti-drag" on one another. 
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Since the liquid velocity vanishes on the plates (as- 
suming no-slip boundary conditions), one might have ex- 
pected suppression of locally induced flows and a cutoff 
in the hydrodynamic interaction at distances much larger 
than w. Such screening is found in q-lD suspensions |^] 
and in the electrostatic analogue of charges located be- 
tween two conducting plates. However, all four diffusion 
coefficients are found to decay as for large p (Fig.lH)). 
This decay is faster than the 1/p dependence in uncon- 
fined suspensions but is still long-ranged — correlations at 
distances ten times w are measurable. Moreover, the ob- 
served decay is slower than the dependence of the 
interaction parallel to a single wall ^ . 

Other particles obstruct the flow induced by the mon- 
itored pair, thus affecting their measured interaction, as 
observed at short distances (Fig. (HI) 14]. One would have 
expected that at large distances the combined contribu- 
tions from numerous particles would shift the pair in- 
teraction by a factor proportional to (j). This correction 
should have been signiflcant already at small (j) because 
of the long-range interaction. However, as can be seen 
in the insets of Fig. O within our experimental accuracy 
there is no concentration effect even at large (j) p^ . 

Since the Reynolds number is very low (~ 10~^), the 
hydrodynamic interactions are carried by linear Stokes 
flows 16j. Together with no-slip boundary conditions on 
the plates and particles the mathematical problem is well 
posed, yet very complicated to solve. The observations 
presented above can be accounted for, nonetheless, on 
simpler grounds. We focus on the flow due to the motion 
of a single particle parallel to the plates. The motion 
perturbs the liquid momentum distribution and displaces 
liquid mass. Far away from the particle one can consider 
only the lowest-order moments of the perturbation — a 
momentum monopole and a mass dipole. (Since mass 
is merely displaced, there is no mass monopole.) The 
momentum monopole undergoes diffusion with absorbing 
boundary conditions at the plates, and its contribution 
to the far fleld is therefore exponentially small in r/w. 
However, the mass-dipole perturbation is not absorbed 
by the boundaries and must propagate laterally. We thus 
expect the far fleld to be that of a two-dimensional mass 
dipole (source doublet), which decays only as The 
entrainment of a second particle by this flow determines 
the pair hydrodynamic interaction, hence the decay 
in the measured correlations. Dipolar flow flelds were 
previously encountered in various q-2D and 2D problems 
|l7l lid Il9j . The above argument clarifles the generality 
of this result. 

Consider a particle located on the midplane and exert- 
ing a force /ix on the liquid. The y axis is taken as the 
second direction parallel to the plates, and z perpendicu- 
lar to the plates [Fig.^b)]. The far velocity fleld, arising 
from a two-dimensional mass dipole oriented along x, is 



given by Va{r) = [fi/ {67rr]w)]Aair), with 
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while Az is exponentially small in r/w. The flow has 
a perpendicular proflle h{z) having a certain value at 
the midplane, h{0) = w^X^ and vanishing on the two 
plates, h{z = ±k;/2) = 0. (In the limit a/w <^ 1 h{z) 
is parabolic and A = 9/16 17].) The dipolar flow pat- 
tern is depicted in Fig. ^fb), where it is contrasted with 
the flow in an unconflned medium [Fig. ^a)]. Due to 
the circulation currents a second particle, located trans- 
verse to the particle motion, will evidently be entrained 
in the opposite direction, hence the negative transverse 
correlation. 



FIG. 4: Flow velocity field away from a particle moving to the 
right (indicated by circle and arrow at the center), in (a) un- 
conflned liquid 16] and (b) liquid confined between two plates 
[Eq. (0)]. Shown are cross-sections of the three-dimensional 
fields along the plane of particle motion. Note the dipolar 
shape of the confined fiow. 



To leading order in a/r, the longitudinal and transverse 
pair interactions are given, respectively, by Ax{rx.) and 
A^(ry). Thus, from Eq. 0, 

A±(p»l) = ±A/p^ A±(p»1) = tA/p^ (3) 

This result is in good agreement with the data of Fig. 
fusing X{a/w :^ 0.45) ^ 0.31. In Fig. |3Jb) we present 
also the results of the full solution in the limit a <^ w 
113, multiplied by 0.31/(9/16) = 0.55 to match the long- 
distance behavior. This calculation reproduces the non- 
monotone behavior of at short distances ^20]. 

We now outline the analysis of the three-body ef- 
fect Particle 1 exerts a force /ix, yet now in the 
presence of particle 2 at r and particle 3 at r', where 
both r, w [Fig. ^b)]. The leading effect of par- 
ticle 3 on the pair interaction between 1 and 2 arises 
from variations of the induced flow over the volume it 
occupies. Because of symmetries of the midplane and 
Eq. (O we flnd that the only derivative contributing 
to this effect is a'^d'^Va/dz'^\z=o ^ {a'^/w'^)va |23. We 
can therefore write the force exerted by particle 3 as 
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fs^a = (67r7^a^/i^^)Ca/3'y/3(r'), where Caf3 are dimension- 
less coefficients wfiicfi, in principle, could be obtained 
from a detailed calculation. From rotational symmetry 
Cxx Cyy Ci^ and C^y -Cyx Ct. The force 
fs perturbs the velocity of particle 2 in the x direction 
by 5vx{v) = {67Tr]w)~^ Aa{T — T')fs^a^ and the perturba- 
tion should be averaged over all third-particle positions 
r^ The inter-particle distances being much larger than 
a, we can assume a uniform probability per unit area, 
^/(7ra^), to find a third particle at r^ Combining this 
with the results for fs and 5vx and rescaling, we find the 
following corrections to the pair interactions: 

5Al^{p) = ±{a/7rw)4>CL J d^p'[A,{p - p')A:,{p') + 
Ay{p - p')Ay{p% (4) 

where ^A^ and SA^ are obtained by setting, respec- 
tively, p = pic and p = py. The integral in Eq. 
however, vanishes. The leading corrections to this result 
that are linear in ^ are found to decay faster than 
[2II and, hence, the vanishing of the three-body term is 
asymptotically exact. 

This approach can be extended to the four-body ef- 
fect j21|J, resulting in finite, small corrections. We find, 
e.g., SA^/A^ = {l/87r)X^{Cl ^ C^){a/w)^\ At our 
largest this amounts to ~ (3 x 10~^)(Cl + C^) which 
(unless Cl,t are unexpectedly large) would be smaller 
than the experimental noise. For most practical purposes 
the hydrodynamic interaction at large distances is thus a 
long-ranged, yet purely pairwise effect even at high con- 
centrations. This result is unique to the q-2D geometry. 

Perpendicular particle motion is restricted in our sys- 
tem by the large a/w ratio and possible repulsion from 
the plates 0|. The flow induced by such motion de- 
cays exponentially with r/w. Increasing the spacing w 
is expected, therefore, to shift the validity range of our 
results to larger inter-particle distances. This effect, how- 
ever, requires further experimental study. Since diffusion 
coefficients are proportional to mobility coefficients, the 
results presented here for Brownian motion should apply 
as well to overdamped driven motion, as in sedimentation 
[2^ or microfluidic transport 2]. 
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